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The production of vector bosons in association with jets contains at least two unrelated scales. 
The first is the mass of the vector boson my and the second is the hard interaction scale giving rise to 
large transverse momenta of the produced jets. The production cross sections depend logarithmically 
on the ratio of these scales, which can lead to a poor convergence in fixed order perturbation theory. 
We illustrate how to resum all leading logarithmic terms using effective theory methods, and show 
that they can be resummed by a simple choice of the factorization scale. Implementing this scale 
choice we show that the large discrepancies between next-to-leading calculations and leading order 
calculations using more traditional choices of scales disappear. 



The production of W and Z bosons in association with 
jets provides one of the most important backgrounds to 
many searches of physics beyond the standard model 
(BSM). This is due to the fact that expected models of 
BSM physics give rise to events containing missing en- 
ergy (from a potential dark matter candidate), jets (from 
cascade decays of new strongly interacting particles) and 
possibly leptons. Since the W and Z bosons can decay 
to neutrinos and/or leptons, they give rise to events with 
the same signatures as the BSM signal one is searching 
for. To distinguish between the expected new physics 
and the Standard Model (SM) background one has to 
study the different kinematical structures of the events. 
For example, the background preferably produces jets 
with small transverse momentum relative to the beam 
axis, while the cascade decays tend to give rise to much 
harder jets. However, in isolating potential BSM signa- 
ture events one often imposes strict kinematic constraints 
to which the background is also subject to. Thus, a de- 
tailed understanding of the differential background spec- 
tra is required in order to successfully separate possible 
BSM signatures from the SM background. 

Much work has gone into detailed calculations of dif- 
ferential distributions for V + jets, where V denotes ei- 
ther a W or a Z boson. Leading order (LO) calcula- 
tions P, [2, H, 0, [E El are available for up to 6 jets in 
the final state, while full next-to-leading order (NLO) 
calculations are only available for one or two jets in the 
final state @, 0, H, El • Recently, first NLO calculations 
for V + 3 jets have become available. In [ll[ the dom- 
inant partonic channels have been calculated at leading 
order in 1/Nc, and the full calculation at leading or- 
der in 1/Nc has followed shortly afterwards The 
subdominant 1/Nc corrections have been argued to be 
numerically small, and first results including all 1/Nc ef- 
fects have been presented recently [13J . Depending on the 
choice of scales, the total cross sections at NLO can differ 
by up to a factor of two from the LO results, while dif- 
ferential spectra can vary even more strongly, especially 
in the important regime of high px jets. 

It is well known that perturbative corrections gener- 
ically contain logarithmic terms that depend on ratios 
of scales in the problem. For the inclusive production 
of a vector boson, the only relevant scale is the mass of 
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FIG. 1: pt distribution for the jet in the pp — > W + j 
process, both at leading order (dashed) and next-to-leading 
order (solid). 



the vector boson my. However, the kinematics of the 
process (and cuts on these kinematic variables) can give 
rise to scales that differ significantly from the mass of 
the produced vector boson. For example, for V + j fi- 
nal states, the transverse momentum px of the jet can 
be far larger than my, and large logarithmic terms of 
the ratio px /my can potentially spoil the convergence of 
fixed order perturbation theory. A simple way to illus- 
trate the size of these logarithmic terms is to compare 
the pt spectrum at LO and NLO. Since the NLO calcu- 
lations contain the logarithms of the form a s log 2 px /my 
and a s \ogpx/my, one expects the difference between LO 
and NLO calculations to get worse as the ratio px/my 
increases. In Fig. [1] we show the px distribution of the 
jet in the LHC environment both at LO and NLO, ob- 
tained using the program MCFM [HG3 

using the same 

cuts as in Ref. [lj] . One can clearly see that the discrep- 
ancy grows with increasing px ■ Note that for V + j final 
states the transverse momentum of the jet is equal to the 
transverse momentum of V, and the ratio of scales can 
therefore be constructed out of the kinematics of the V 
boson alone. 

For processes with at least two jets in association with 
the produced vector boson, the individual transverse mo- 
menta of the jets are no longer related to the transverse 



2 



(WWV 



fWWV 



* 000000 ■> '000000 4 



FIG. 2: Feynman diagrams contributing to the process pp 
W++j. 



AA/WV 



= C 



a 



FIG. 3: Matching of the amplitude for the partonic subprocess 
qq — > Vg. A similar relation holds for qg — > 



momentum of the vector boson. Since the transverse mo- 
menta of the jets can balance each other, the pt of each 
jet can get large independent of the transverse momen- 
tum of the vector boson. However, the same arguments 
go through as for V + j final states, and logarithmic de- 
pendence on the ratio pr / my is expected again. 

At higher orders in perturbation theory more and more 
powers of the large logarithms can arise. In general, for 
any ratio of scales r one finds up to two powers of log- 
arithms for each power of a s , and for sufficiently large 
ratios of scales r at least the leading logarithmic (LL) 
terms of the form a™ log 2 ™ r have to be resummed to all 
orders to obtain a reliable prediction. 

Effective theory methods naturally separate widely dif- 
ferent energy scales from one another, and they are ide- 
ally suited to understand the corresponding logarithmic 
terms. Effective field theories are designed to reproduce 
an underlying theory at long distances, while short dis- 
tance physics are captured in Wilson coefficients, which 
have to be adjusted order by order in perturbation theory 
to match the full theory at short distances. 

The power of effective field theory techniques is that 
for processes that involve several widely separated energy 
scales Ai ^> A 2 ^ . . . A„ these scales can be removed 
one by one by matching onto successive effective theo- 
ries, where each effective theory is valid below one of the 
scales Ai in the problem. The Wilson coefficients in gen- 
eral depend on a renormalization scale and logarithms 
of the ratio = p,i/Ai will appear in the expressions. 
However, since each matching calculation only depends 
on the single scale A^ (larger scales have been removed in 
a previous matching step, and the matching calculation 
is independent of the smaller scales) , the choice Hi ~ Aj 
minimizes these logarithmic terms. One can then use the 
renormalization group equations (RGEs) to sum the log- 
arithms of ratios fii+x/fii between the different matching 
scales. 

Since processes observed at high energy colliders typ- 
ically contain light particles with large energies, the ap- 
propriate effective field theory which reproduces all long 
distance physics of the standard model (SM) is soft- 
collinear effective theory (SCET) 0, OJ, EE \lk ■ SCET 
contains collinear fields to describe the light, energetic 
particles, which can only interact with one another via 
the exchange of soft particles. For details on the con- 
struction and applications of SCET we refer the reader 
to the literature. 

We start by considering the process pp —* V + j , for 
which the tree-level amplitude can be calculated from 
the diagrams shown in Fig. O We try to separate the 



two scales Ai = px and A 2 = my by using successive 
matching calculations. We will keep our discussion of 
the matching calculation very schematic, and ignore all 
complications, for example those due to operator mixing. 
This is sufficient to understand the resummation of the 
leading logarithmic terms at tree level, which is the aim 
of this letter. However, the procedure can easily be gen- 
eralized in a completely straightforward manner to allow 
the inclusion of higher order corrections. 

The intermediate propagator in Fig. [2] scale as l/p T 
for large values of p\, and only gives rise to short dis- 
tance physics. Thus, the same amplitude is reproduced 
in the effective theory using operators containing only 
the incoming partons and the outgoing parton and vec- 
tor boson, and the effect of the intermediate propagator 
is contained in the Wilson coefficients Q of these op- 
erators, as is illustrated in Fig. [3] Note that since we 
consider events with px 3> my, we can treat the vector 
boson as massless and as part of SCET. As argued above, 
the logarithms in the matching coefficients are minimized 
if one chooses fix ~ Ai = pt in this calculation. Thus, 
we find for the amplitude 

A ~ C 0qq (p T )(O qq ) PT + C 0qg (PT)(O qg ) PT , (1) 

where {O qq ) denotes the matrix element of the operator 
shown on the RHS of Fig. [3J and O qg denotes a similar 
operator for qg initial states. Note that the matrix el- 
ements of the operators Oij are evaluated at the scale 
fi = px, as indicated by the subscript. The matrix ele- 
ment of the operator at a different scale can be obtained 
using the RGE of the effective theory, 

(O ij ) PT =U 0tj (pT,lj)(O ij ) li , (2) 

and the evolution kernels C/q can be calculated straight- 
forwardly in perturbation theory. 

The next scale to integrate out is the electroweak scale, 
and with it my . Since below this scale we can no longer 
treat the vector boson V as a propagating particle in the 
effective theory, we have to remove V from the theory 
altogether. This can be accomplished by considering the 
forward scattering matrix element as shown in Fig. 3] and 
integrating the intermediate V boson out of the theory. 
The corresponding Wilson coefficients Cj? t will depend 
on my, and the associated logarithms can be minimized 
by choosing /i 2 ~ A 2 = my. Thus, the differential cross 
section is given by 

dcr ~ C O n {VT)V 2 o qq (pt, m v )Cjr qq {m v ){T qq ) mv 

+ c o qt ,(PT)Uo qq (PT,my)Cjr iig (my)(F qg ) mv , (3) 
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FIG. 4: Matching of the forward scattering matrix element 
for the partonic subprocess qq — > Vg. A similar relation holds 
again for qg — > Vq 



where (fqq) denotes the matrix element of the operator 
shown on the RHS of Fig. [H and (J- qg ) is the obvious 
generalization to qg initial states. 

Finally, we have to calculate the matrix elements of 
the remaining operators. This can be accomplished us- 
ing the factorization properties of QCD and SCET, such 
that the remaining long distance physics is absorbed into 
the parton distribution functions of the incoming partons 
and the jet functions of the outgoing partons. For suffi- 
ciently inclusive jet definitions, the jet properties can be 
calculated entirely in perturbation theory, and one writes 



U Fu {m v ,li)Eij(ji)fi{jx)Sj(ji) ■ (4) 



The scale /x denotes the scale at which the parton distri- 
bution functions are evaluated, the Kernel U^Amy, A*) 
denotes the RG evolution of the operator from the scale 
my to the scale \x, and Eij denotes extra perturbative 
terms that depend on the scheme used for the parton 
distribution functions. 

Combing these results, and using the fact that by con- 
struction the combination of all matching coefficients re- 
produces the full partonic cross section expressions, the 
result reads 



da ~ ^ ( X i > X j)fi ( X i> M)/j ( X j > A*) 

ij 

x u o t . (pr, m v )U Ti . (my, . 



(5) 



All logarithmic terms of the ratios my/px and [if my 
have been resummed in the kernels Uo i} (pr> uiy) and 
CT^r (my, /it). However, the result in ([5]) can be simplified 
further by using a simple relation between the two evolu- 
tion kernels, which is correct at LL accuracy. The leading 
logarithmic structure is governed by the so-called cusp 
anomalous dimension of the operators Oij and Tij , which 
can be calculated in SCET from collinear one-loop dia- 
grams involving only strongly interacting collinear fields 
in a given direction. Since the operator Tij is obtained 
from the forward matrix element of two operators Oij 
by integrating out the vector boson which is not strongly 
interacting, the evolution kernels are related at LL order 
by the simple relation 




FIG. 5: Sample Feynman diagrams for initial-state radiation 
and final-state radiation contributions to pp — > V + 2 jets. 



Since furthermore the evolution kernels satisfy the simple 
relation U(fii, H2)U(fi2, A*a) = U(/j,x, /is), we can simplify 
our result as 
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Thus, by choosing to evaluate the parton distribution 
functions at a dynamical scale \x = pt, the dependence 
on the evolution kernel disappears entirely, and for that 
scale choice the result is just given by the usual fixed 
order expression. 

This result is not unexpected, since it simply states 
that the relevant scale for events with px 3> my is given 
by the large transverse momentum in the process, rather 
than the mass of the vector boson. In order to cover 
processes for which px <C my (complementary to our 
previous assumption) , for which the relevant scale should 
be fj, ~ my, several analyses have used the dynamical 
scale fi = [my +p Tv } 1/2 [3, H, H EH . Since for V + j 
final states the transverse momentum of the vector boson 
Pt v is equal to that of the jet, our result shows that this 
scale setting does indeed sum the large logarithms of the 
form \ogmy / pt ■ 

For processes with two or more jets in the final state, 
large logarithmic terms will originate from phase space 
integrations. However, these logarithms can still be re- 
summed using a very similar analysis as performed for 
pp — > V + j case. To show this we study the kinematics 
for pp — > V + 2j. Many Feynman diagrams contribute to 
this process, but for large partonic center of mass energy 
the dominant contribution to the cross section arises from 
diagrams where the V boson is radiated from an initial or 
final quark with either soft momentum or collinear to the 
quark it is radiated from. This gives rise to effects that 
are enhanced by a factor of px/my. Two representative 
Feynman diagrams of such initial or final state radiation 
are shown in Fig. [5] Both of these processes result in 
kinematics where the two jets are back-to-back in the 
transverse plane, with the angle 4>jj between them near 
7T. We illustrate this effect in Fig. [SJ where the distribu- 
tion in 4>jj for different values of a cut on the smallest 
available pr of the jets is shown. One can see that the 
larger the cut on pt, the more the distributions peak to- 
wards 4>jj — 7r, in agreement with our expectations. Note 
that there is another kinematical region that is enhanced 
by a large propagator, which is when the two jets are 
collinear to each other, and are recoiling against the V 
boson. It is the size of the jets which regulates this singu- 
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FIG. 6: Distribution in (fijj for pp — > W + 1j and for different 
cuts on the minimal jet transverse momenta. 

larity, and in fact the small feature near 4>jj = 0.4 is due 
to this kinematical configuration (AR = 0.4 was used in 
the jet algorithm). However, this region is enhanced by 
only a power of 1/(AR), and for large enough values of 
px the initial and final state radiation of the V boson 
dominates. 

To resum the logarithms of order my/px, we again 
perform a two step matching calculation. In the first 
step, we integrate out the propagators with virtuality of 
order p Tl resulting in a matching calculation as indicated 
in Fig. [7J In order to avoid large logarithmic terms in this 
matching calculation, we choose n ~ px- Note that for 
the dominant initial and final state radiation contribu- 
tions the V boson does not participate in this matching 
calculation. In a second matching step, the V boson, 
together with all other propagators with virtuality of or- 
der my are integrated out. The natural scale for this 
matching is again /i ~ my. Finally, the matrix elements 
are calculated, resulting in parton distribution functions 
evaluated at that final scale \x. In order to connect the 
different scales to one another, RG evolution is used, and 
using a similar argument as in the case of V + j final 
states, one can show that to LL order the evolution ker- 
nel Uo connecting the scales px to my is related to the 
kernel Up describing the running from my to /i by ([6]). 
Thus, all leading logarithms can again be summed by 
choosing the factorization scale to be of order pt of the 
jets. Note that despite the derivation being so similar as 
for the case of V plus a single jet, to our knowledge the 
scale choice derived in this work has not been used in any 
published results 1 . 

Unlike the change from n = 1 to n = 2 in pp — > V + n 



FIG. 7: First matching step at the high scale /j, ~ pr- 



jets processes, a further increase to n > 3 does not give 
rise to qualitatively new kinematic situations, in that in- 
dividual transverse jet energies can far surpass the vector 
boson's transverse energy. For any given set of outgoing 
4-momenta, there are of course more scales in the prob- 
lem, but for generic events, where each of the jet energies 
scales like their respective px, the very same framework 
based on EFTs can be applied. 

We now study the effect of the resummation of the 
leading logarithmic terms by comparing the LO calcu- 
lation with the new choice of scale (J, ~ px with the LO 
calculation using the traditional scale choice fi ~ my and 
the NLO calculation with fi = my. Note that our scale 
setting procedure only resums the logarithms correctly 
at LO/LL accuracy, and we therefore do not change the 
scale used in the NLO calculation. To be concrete, we 
chose 2 /x 2 yn = (m.hadr/2) 2 + m\ , where m hadr is the dijet 
invariant mass of the hadronic system. The scale depen- 
dence at NLO is considerably reduced over the depen- 
dence at LO, such that the choice of scale used in the 
NLO calculation does not affect our conclusions. We will 
present results for W + jets in the final state, but we 
have checked that our results apply equally well for Z + 
jets in the final state. 

In Fig. [5] we show on the left the distribution of the 
transverse momentum of the softer jet in pp — > W + + 2j 
events. For small values of px the two different choices 
for the scale give comparable results, with both results 
agreeing reasonably well with the NLO calculation. For 
larger values of px the LO calculation with the traditional 
choice of scale setting fi — my disagrees more and more 
with the NLO result, while the LO calculation with the 
dynamical scale setting agrees with the NLO results over 
the full range of px- To illustrate that this is not just 
an artifact of the distribution in transverse momenta, we 
also show the Hx distribution, where we have imposed a 
cut px > 100 GeV to select events where the difference 
between the two scale setting choices becomes more pro- 
nounced. Again, we see that that the LO calculation for 
the traditional choice of scale setting /i = my does not 
agree well with the NLO calculation, but that the new 
scale choice improves the situation significantly. We have 
compared the traditional and new scale setting for many 
other observables, and in all cases the new choice of scale 



1 An observations that the traditional scale choice in LO calcula- 
tions gives poor convergence with NLO calculations has recently 
been made |22| , together with a suggestion for using a dynamical 
scale similar to ours in spirit. For vector boson fusion, dynamical 
scales have also been found empirically to improve the conver- 
gence of perturbation theory 1231 . 



2 This scale choice is somewhat ad hoc, but consistent with its 
scaling with px- Other choices which satisfy the correct scaling 
with px differ at most by logarithms which are not enhanced in 
the limit my/px — -> 0. 
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FIG. 8: Comparison between the two scale settings fi — my (dashed, blue) and /i = /idyn (dotted, red) and the NLO prediction 
(solid, black) for the process pp — > W + + 2j. Left: pr spectrum of the second jet with a cut of pr > 20 GeV. Right: Ht 
spectrum with pr > 100 GeV . For high-pr events the NLO calculation becomes less reliable as statistical fluctuations gain 
more prominence. 



improves the distributions. 

In Figure [9] we compare for n — 3 the LO predictions 
for the pT(j,min) spectrum - the smallest of the three 
transverse jet momenta - with traditional scale setting 
and with the dynamical scale setting. We again observe 
that the net effect of resummation is a lowering of the 
differential cross section, and that the region of large pr 
receives a larger relative correction than the small pt 
region, i.e. a correction of the shape of differential cross 
section. This agrees qualitatively with the heroic fixed- 
order NLO calculation recently performed [HI, [HJ , both 
for the pt and similarly for the Ht distributions. 

It should be noted that the simple scale setting proce- 
dure used above gives results at leading order in pertur- 
bation theory with the leading logs resummed. The ef- 
fective field theory methods introduced, however, can be 
used to improve these results by including higher orders 
in fixed order perturbation theory and/or higher order 
terms in the resummation calculation. NLO corrections, 
for example, can be included by calculating the matching 
coefficients Co id , C.Fy and to higher orders in per- 
turbation theory, while higher orders in the resummation 
can be included by computing the evolution kernels Uo i:j 
and Ujr. . to higher orders. Note that we expect the rela- 
tion between the two evolution kernels given in (|6|) to be 
violated at NLL order, such that the NLL resummation 
can not be accomplished by a simple choice of scales in 
the parton distribution functions. 



In this letter we have investigated the effect of leading- 
logarithmic resummation for jet production with associ- 
ated vector boson (V = W, Z) production at the LHC. 
The resummation is achieved via a series of matching 
onto versions of soft-collinear effective theory and run- 
ning of the coefficient functions to their natural scales 
using their rcnormalization-group properties. The result 
simplifies to the known partonic cross section expression 
folded with the appropriate parton distribution functions 
if the parton distribution functions are evaluated at a 
high scale which scales like the transverse momentum of 
the jets. The main purpose of resummation is to improve 
the convergence of the perturbative expansion which we 
have demonstrated specifically for events with high trans- 
verse momenta in pp — * V + n jets processes. 
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